Friction coefficient for deep-inelastic heavy-ion collisions. 



G.G.Adamian 1 ' 2 , R.V.Jolos 1 and A.K. Nasirov 1 ' 2 
1 Joint Institute for Nuclear Research, Dubna, 14-1980 Russia 

A.I.Muminov 

2 Heavy Ion Physics Department, Institute of Nuclear Physics 702132 Ulugbek, Tashkent, 

Uzbekistan 
(February 9, 2008) 

Abstract 



Based on the microscopic model, the friction coefficient for the relative 

motion of nuclei in deep- inelastic heavy-ion collisions is calculated. The radial 

dependence of the friction coefficient is studied and the results are compared 

with those found by other methods. Based on this result, it was demonstrated 

that the kinetic energy dissipation in deep-inelastic heavy-ion collisions is a 

gradual process which takes up a significant part of a reaction time. An 

advantage of the suggested method is that it allows one to consider the relative 

motion of nuclei and the intrinsic motion self-consistently. 
25.70.Gh, 25.70.Jj 
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I. INTRODUCTION 



Nuclear friction is an important ingredient of theoretical approaches to a variety of 
nuclear physics phenomena, such as dynamic thresholds for compound nucleus formation 
[ jH|, enhancement of neutron emission prior to fission 0,0], width of mass and charge 
distributions in deep-inelastic heavy-ion reactions M, and the width of giant resonances 
||. There are many experimental results on deep- inelastic heavy- ion collisions (DIC) and 
fusion-fission reactions which need the introduction of the nuclear friction concept for their 
interpretation. This stresses the importance of understanding the nature of nuclear friction. 

The present paper is devoted to calculations of the friction coefficient for DIC heavy-ion 
collisions. Its appearance is stimulated not only by the possibility to perform more exact 
calculations than earlier, but also by the new experimental results which require a more 
detailed microscopic theory for their interpretation than was necessary before. 

Different theoretical approaches to this problem are known. The majority of them are 



based on the assumption that the dissipative mechanism is of a one-body nature [10-12 



These models differ in the structure of the intrinsic excitations that are taken into account. 



These can be pure collective excited states ||T3|-[T6| or incoherent particle-hole excitations 
T0| , p!8|JT9| . Some of the models do not include nucleon exchange and consider only particle- 



hole excitations, with both the particle and the hole belonging to the same nucleus or 
[Toll . The models also differ in the approximations they use for including the finite 



vice versa 



decay time of one particle-one hole (lp-lh) excitations in more complicated configurations 
(2p-2h, and so on). Many approaches 0.[[|.[l^.|l[] implicitly use the statistical assumption 
of rapid equilibration of the noncollective intrinsic degrees of freedom and, therefore are not 
applicable to the description of the initial phase of the reactions, where the main part of 
kinetic energy dissipation takes place. 

The contribution of the actual lp-lh state or more complicated ones to the dissipation 
process depends on the occupation numbers of the single-particle states and their evolution 
during the reaction. However, in the calculations of the nuclear friction coefficient performed 
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up to now, the statistical assumption on the excitation energy distribution was realized 
in the usual way, meaning the introduction of temperature and, correspondingly, of the 
Fermi occupation numbers at the very beginning of the reaction. In principle, temperature 
introduced in this way is a time-dependent quantity. However, in practice temperature was 
considered to be a constant corresponding to the total excitation energy ]To]Ji7| . Thus, the 
time dependence of the single-particle occupation numbers was not taken into account. Only 



in the approach based on the Dissipative Diabatic Dynamics (DDD) |22|j23|| was the evolution 
of the single-particle occupation numbers taken into account, but under the assumption of 
diabaticity. There are some doubts however, about the validity of the DDD concept. It 
is also known from the calculations of inelastic processes in nucleus-nucleus collisions that 
appreciable energy dissipation takes place even before the first crossing of the single-particle 
levels near the Fermi surface |24[ . 



We should also mention the approach suggested in ||25|| , where relative and intrinsic 



motion were consistently treated in a time-dependent theory of heavy-ion collisions. The 
authors presumed neither weak coupling between the relative motion and the intrinsic exci- 
tations nor the canonical distribution function for the density operator of intrinsic motion. 
However, the analytical expressions for the friction tensor and other characteristics of the 
energy transport obtained in |25| were not applied to calculate them. 



Thus, it is the aim of the present paper to take into account the time evolution of 
the single-particle occupation numbers during the reaction by numerical solution of the 
master equation for them and, based on this result, to perform calculations for the friction 
coefficient. Since the occupation numbers found in this way correspond to the current 
kinetic energy losses, this means that the relative and intrinsic motions are considered self- 
consistently. Our model makes it possible to take into account explicitly the influence of the 
nuclear shell structure on the collision process. Moreover, we improved the single-particle 
approximation by a phenomenological allowance for the residual interaction which is treated 
in the so called r-approximation. The radial friction coefficient is calculated as a function 
of the mass and charge of the reaction participants. 



The general formalism is given in Section 2. The results of the calculations are presented 
in Section 3. A summary is given in Section 4. 



II. BASIC FORMALISM 

It is convenient to start with the total Hamiltonian of a dinuclear system written in the 
form 

H = H rel (R; P) + H in (0 + 5V(R, £), (1) 

where the Hamiltonian of a relative motion, 

ff rd (R;P) = ^+V(R), (2) 

consists of the kinetic energy operator and the nucleus-nucleus interaction potential V(R). 
Here, R is the relative distance between the centers of mass of the fragments, P is the 
conjugate momentum, and fi is the reduced mass of the system; £ is a set of relevant 
intrinsic variables. The last two terms in (|IJ) describe the internal motion of nuclei and the 
coupling between the relative and internal motions (for details, see 1 26,271). It is clear that 



the coupling term leads to a dissipation of the kinetic energy into the energy of internal 
nucleon motion. 

Neglecting at the moment the residual nucleon-nucleon interaction, whose effect will be 
included later, we take a sum of the last two terms (|l|) as a single-particle Hamiltonian of a 
dinuclear system 

H in (0 + 5V(R,0 = H(B.(t),£) + h residuah 

= E + v P (n - *(*)) + M*)) , (3) 

where m is the nucleon mass and A = Ap + Aj- is the total number of nucleons in the system. 

Then, in the second quantization representation, the Hamiltonian 7Y(R(t),£) can be 
written as 
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W(R(t),0 = $> P a£a P + 5> T a T a T + ^(R^a+a^, (4) 

P T iji%> 

where 

V r ^(R(t))a+a i , = 51 App/(R(t))a£a P , + ^ A^(R(*))ara r , + ( 5 ) 

i^i' P^P' T^T' 

£g PT (R(t))(a+a r + h.c.). 

T,P 

Here P = (n P ,j P , l P , m P ) and T = (n T ,j T , l T , m T ) are the sets of quantum numbers charac- 
terizing the single-particle state in an isolated projectile and the target nuclei, respectively. 
The single-particle basis is constructed from the asymptotic wave vectors of the single- 
particle states of the noninteracting nuclei-the projectile ion \P > and the target nucleus 
|T > in the form 

\P> = \P > -- Y \ T >< T \ P >> ( 6a ) 

\f> = \T > -- Y \P >< P\T > ■ (6b) 
2 p 

For this basis set, the orthogonality condition is satisfied up to terms linear in < P\T >. 
Then 

AgL(R(t)) =< P\V T {v)\P' >, (7a) 

A©(R(t)) =< T\V P (r - R(t))|T' >, (7b) 

gprW)) = \< P\Vp(r - R(t)) + V T (r)\T > . (7c) 

The nondiagonal matrix elements Ap z P ,(A^,) generate the particle-hole transitions in the 
projectile (target) nucleus. The matrix elements gpT are responsible for the nucleon exchange 
between reaction partners. These matrix elements were calculated using the approach pro- 
posed in p8| , |29| . In (|j), e P {r) are the single-particle energies of nonperturbed states in the 
projectile (target) nucleus. The coupling between the intrinsic nuclear degrees of freedom 
and the collective variable R is introduced by the R dependence of the sum of the single- 
particle potentials in ([3]). Since the trajectory calculation shows that the relative distance 
R(t) between the centers of the interacting nuclei could not be less than the sum of their 
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radii the tail of the partner single-particle potentials can be considered as a perturbation 
disturbing the asymptotic single-particle wave functions and their energies. 

It is convenient to include the diagonal matrix elements of Vai(H(t)) in Hi n , introducing 
the renormalized R(t)-dependent single-particle energies 

e P (R(t)) = e P + < P\V T (r)\P >, (8a) 
e T (R(*)) = e T + < T\V P (r - R{t))\T > . (8b) 

When the nuclear forces begin to act between the colliding nuclei, the velocity of their 
relative motion can be considered as a small quantity compared to the Fermi velocity. Then 
the speed of the nucleons is mainly associated with their intrinsic motion. Since the relative 
(collective) motion is rather slow compared to the intrinsic motion, the perturbation of 
the intrinsic motion produced by changing the coupling to the relative motion (R) can 
be assumed to be small during some small time interval At of an arbitrarily chosen time 



t [13]. The small parameter in our consideration At thus characterizes the time interval 
during which the R-dependent mean field of the combined dinuclear system changes so little 
that we can neglect the effect of this changing on the intrinsic motion. At the same time, 
the characteristic time At can not be taken smaller than the relaxation time of the mean 
field. The situation described above is suitable for applying the linear response theory to 



a description of dissipative heavy- ion collisions [O]. For this reason, we start from the 



expression for the friction coefficient of the radial motion obtained in that approach ||13| , 

Bg\t). (9) 



lR R (R(t)) = £ 

ik 



dV ik (K) 2 



OR 



B%\t) = \\\ dt'^^- exp [tJ\ Bin[a w (R(0)(t - 0] 
n Jt-At nl \ Tiic i 



2 ft jM (t' -t) n _ _ ft' -f 

Tik 

x [n k {t') - hi{t% (10) 



where r^- = TiT^jiji + Tfc); Ti is the parameter describing the damping of the single-particle 
motion. The expression for Tj is derived in the theory of quantum liquids pO| , pT| (see 
Appendix A); faoij = £i(R) — £j(R) is the energy of the single-particle transition in one of 



the nuclei as well as between the interacting nuclei. The important ingredients of this formula 
are the occupation numbers of the single-particle states fii{t). Since the excitation energy of 
the interacting nuclei changes significantly during the course of the collision, it is necessary 
to take into account the time dependence of the occupation numbers. The importance of 
this point was already stressed in [|T^]. At the same time, new experimental data indicate 
that the assumption of the fast statistical equilibration of the excitation energy during the 
collision time, i.e., an introduction of a time-dependent temperature and Fermi occupation 
numbers is not adequate for the physical picture. As already mentioned in the Introduction, 
the calculations of 7^ performed up to now have been done under the assumption that the 
occupation numbers can be taken as the Fermi occupation numbers 



n 



j 



;i + exp((E,-A)/0))- 1 , 



where is the temperature corresponding to the total excitation energy of a dinuclear 
system. 



To find the time- dependent occupation numbers rij(t), we developed in ^,^,0 a 
method which is described briefly below for completeness of the presentation. 

Since explicit allowance for the residual interaction requires extensive calculations, it 
is customary to take the two-particle collision integral into account in linearized form (r- 
approximation) . 

Then, the equation for the single-particle density matrix n takes the form 

iU 3 ^- = [H(R(t))M)) ~ jiHt) ~ n\n(t))] , (11) 

where n eq (H(t)) is a local quasi-equilibrium distribution, i.e. a Fermi distribution with the 
temperature T(t) corresponding to the excitation energy at the internuclear distance R(i). 



Substituting our Hamiltonian (|]) into (11 , we have 



dn (t) 

= E [V ik (R(t))n ki (t) - V ki (R(t))n ik (t)\ 

--[nm-n?{t)} , (12) 

n 



where Hi is a diagonal matrix element of the density matrix. The approximate equation for 
nondiagonal matrix elements takes the form 

2i. 



.+ dn ik (t) 
in — — — = n 



dt 



w«(R(t)) 



Tik 



n%k(t) 



+ V ki (R(t)) [h k (t) - hi{t)} 



(13) 



where we have used the notations = [si — e k ) /h. 

Assuming incoherence in the phases of the nondiagonal matrix elements, we use the 
following approximation to simplify equation ([13|) : 



^2V k ,i{R{t))n k , k {t) -^V^nu'it) w V ki (R(t)) [n k {t) - m(t)) 



As formulated above, we shall consider the solution of equations ([12]) and (|T3| ) for a small 
time interval At of an arbitrarily choosen time t. Then the solution of equation ([13]) can be 
written as 



n ik (t) 



in 



dt'V ik (n(t')) exp < i J dt" 



Q ki (R(t") + 



Tik 



x \n k (t') -hi{t')} 



(14) 



where t < t < t + At. Substituting this result into Eq. (|T^) and transforming this equation 
to an integral, we obtain 



hi{t) =exp ( l — - j | hi(t) + i J dt'hf (R(t')) exp - 



t t' 



/ dt> J dt"Q ik {t',t") exp 



t t 



't" -t" 



Tik 



where 



QikM = ^Re { V ik (R{t))V ki {R{t')) exp 



2 | d*"a) W (R(*")) 



The formal solution of Eq. (|T5| ) is 



n l {t)=n?{n{t)) 



1 — exp 



-At 



Ti 



rii(t) exp 



-At 



Ti 



(15) 
(16) 



(17) 



where 

t 



n,(f) = n t (t) + dt'n ik (t', Sm[ ^ ( R (t ' )} MO " (18) 
In fact, Eqs. ([17]) and dl8|) present an integral equation for rij(i). 



III. RESULTS AND DISCUSSION 

In this section, we present the results of calculating the radial friction coefficient and the 
kinetic energy losses as a function of the inter-nucleus distance for trajectories corresponding 
to DIC. The initial projectile energy, atomic masses, and charges of the colliding nuclei are 
the initial information used in the calculations. The single-particle potentials of the colliding 
nuclei are taken in the Woods-Saxon form with the parameters ro=1.18 fm and a = 0.54 
fm. The characteristic time parameter At is taken to be equal to (0.8-1.0)-10~ 22 s. 

The equations of motion for R(t) and the single-particle occupation numbers fii{t) have 
been solved numerically, transforming differential equations into finite difference equations 
with the time step At and the initial conditions i?(0)=20 fm and fii (0) = 1 or for oc- 
cupied and unoccupied states of the noninteracting nuclei, respectively. Matrix elements 
AgL, (A&l), and gPT have been calculated using the procedure developed in |p6|,p7| . 

The relative motion along the trajectory depends on the nucleus-nucleus interaction 
potential, which is determined by a double folding of the effective nuclear and Coulomb 
interactions of the nucleons with the nuclear densities of the interacting nuclei. Because of 
nucleon exchange and particle-hole excitations, the nuclear densities of the colliding nuclei 
evolve during the reaction and the nucleus-nucleus potential changes correspondingly. This 
effect is included in our calculations. As an example of these calculations we present the 
radial friction coefficient for the 64 Zn(440 MeV) + 196 Pt reaction as a function of R for an 
approach phase of the reaction (solid curve) in Fig. 1. The results of the classical model 
of D.H.E. Gross and Kalinowski [[33| (dashed curves with stars) are shown for comparison, 
along with the results of the microscopic model developed in [Il5j which are obtained with 



S 



a constant temperature (dashed curve without stars). The curves correspond to different 
temperatures which increase from 0.5 MeV (bottom curve) to 2 MeV (upper curve). It can be 



seen that the difference between our results and the results of |33[] increases with the increase 
in the overlap of the colliding nuclei. Comparing our results with the results of we can 
see that our radial friction coefficient coincides with those found in [13] for temperatures 



increasing as R decreases. This kind of behavour of the radial friction coefficient obtained 
in ||17|| . It is natural since it corresponds to an increase in the kinetic energy loss during 
the approach phase of the reaction. Unfortunately, the results for larger values of R are 
not presented in ||15|| . Thus, our calculations demonstrate the importance of inclusion of the 
time dependence of the single-particle occupation numbers in the calculations. 

In Figs. 2 and 3, the radial friction coefficients 7^ are given as a function of R for both 
approach (solid curve) and departure (dashed curve) stages of the 64 Zn(440 MeV) + 196 Pt 
and 56 Fe(480 MeV) + 208 Pb reactions. It can be seen that ^rr takes a larger value for the 
approach stage compared to the departure stage. For the second reaction, this difference is 
significant. 

In Fig. 4 and 5, the dependence of the total kinetic energy losses as calculated by our 
method (solid curve) and as based on the classical model of realized with the code 
TRAJEC H is demonstrated for the 64 Zn(440 MeV) + 196 Pt and 56 Fe(480 MeV) + 208 Pb 
reactions, respectively. The solid and dashed arrows indicate the moment corresponding to 
the turning point of the trajectory in our calculations and in the classical model, respectively. 
It can be seen that, in contrast to the results of the classical model, where the kinetic 
energy is dissipated during a short time interval of the order of 0.4 -10 -21 s at the very 
beginning of the reaction, this process takes significantly larger time in our model. This 
result is a consequence of the smaller value of our friction coefficient relative to those used 
in the classical phenomeno logical models. As seen in Fig.l, the friction coefficient obtained 
in [|3i| significantly exceeds our friction coefficient where the interacting nuclei strongly 
overlap. Nevertheless, our friction coefficient reproduces the total amount of kinetic energy 
loss observed experimentally. If the results of the classical model calculations support the 

9 



idea of the fast kinetic energy loss and thermalization, our calculations support the idea of 



gradual kinetic energy dissipation as was argued in [||,|3^]. 

In Fig. 6, we show a correction to the nucleus- nucleus interaction potential connected 
with the rearrangement of nuclear densities during the approach phase of the reaction for the 
collision of 64 Ni(320 MeV) + 208 Pb. It can be seen that the correction increases in absolute 
value with a decrease of the internucleus distance R. For the same reaction, the nucleus- 
nucleus interaction potential found in the sudden approximation (dotted line) is shown in 
Fig. 7 with the dynamical correction (dashed line), which is discussed just above, and the 
kinetic energy as a function of the trajectory (solid line). We show that the correction to 
the sudden approximation produced by a rearrangement of the particle distribution in the 
interacting nuclei is important and must be included in the consideration. 



IV. SUMMARY 

In conclusion, we have calculated the friction coefficient for DIC based on the microscopic 
model of the structure of the colliding nuclei and avoiding the assumption of a fast statisti- 
cal equilibrium of the dissipated kinetic energy. Our results demonstrate the importance of 
considering the friction coefficient as an exact dynamic function of the single-particle occu- 
pation numbers without the assumption of the fast statistical equilibrium of the excitation 
energy. Based on this result, we have demonstrated that the kinetic energy dissipation in 
DIC is a gradual process which takes up a significant part of a reaction time. 
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APPENDIX A: 



The value of r, is calculated using the results of the theory of quantum liquids p0|j31 



V2' 



71 



T, 



(a) 



(fK-g) 2 + -(fK + gf 

» 



(tt&k) 2 + (ii-X 



(") 

K 



(Al) 



where 



Q K (t) = 3.46, 



Em 



\ < A K (t) > 

is the effective temperature determined by the amount of intrinsic excitation energy E* K = 
E*j} z) + E^ N) - < A K (t) >=< Z K (t) > + < N K (t) >, \P(t), and E*^ a) (t) are the mass 
number, chemical potential, and intrinsic excitation energies for the proton (a = Z) and 
neutron (a = N) subsystems of the nucleus K(K = P,T), respectively, (for details, see 
I. Furthermore, the finite size of nuclei and the available difference between the numbers 



of neutrons and protons need to use the following expressions for the Fermi energies |]3T 

2 



e F K - £ F 



■F K 



e F 



[l + 2f 



1 + 3(1+2/ 



< N K > - < Z K > 

< A K > 
n < N K > - < Z K > 

< A K > 



(A2) 



where €f=37 MeV, 



f K fin 
fx = fin 



< A K >V3 



2/3 {fin fex)i 



< A K >V3 



if ~ f ) 

\J in J ex/ 



(A3) 



and /j n =0.09, //„=0.42, / e:r =-2.59, /g X =0.54, (?=0.7 are the constants of the effective nucleon- 
nucleon interaction. 
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FIGURES 

FIG. 1. The radial friction coefficient calculated according to the present model for the approach 
phase of the Zn(440 MeV) + 196 Pt reaction as a function of R (solid curve). For comparison 
the results of the classical model of D.H.E. Gross and H. Kalinowski [32] are shown (dashed curve 
with stars) and also the results of the microscopic model developed in [15] which are obtained 
with a constant temperature (dashed curves without stars). The curves correspond to different 
temperatures which increase from 0.5 MeV (bottom curve) to 2 MeV (upper curve). 

FIG. 2. The radial friction coefficients jrr as a function of R for both the approach (solid 
curve) and departure (dashed curve) stages of the 64 Zn(440 MeV) + 196 Pt reaction. 

FIG. 3. The same as in Fig.2 but for the 56 Fe(480 MeV) + 208 Pb reaction. 

FIG. 4. Dependence of the total kinetic energy losses calculated by our method (solid curve) 
and based on the classical model of [32] realized with the code TRAJEC [33] for the 64 Zn(440 
MeV) + 196 Pt reaction. The solid and dashed arrows indicate the moment corresponding to the 
turning point of the trajectory in the present calculations and in the classical model, respectively. 

FIG. 5. The same as in Fig.4 but for the 56 Fe(480 MeV) + 208 Pb reaction. 

FIG. 6. Correction to the nucleus-nucleus interaction potential connected with a rearrangement 
of the nuclear densities during the approach phase of the reaction for the collision of 64 Ni(320 MeV) 
+ 208 Pb. 

FIG. 7. The nucleus-nucleus interaction potential found in the sudden approximation (dotted 
line) and with the dynamic correction (dashed line) presented in Fig. 5, and the kinetic energy as 
a function of the trajectory (solid line) for the 64 Ni(320 MeV) + 208 Pb reaction. 
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